Present and Future Gamma-Ray Probes of the Cygnus OB2 Environment 
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The MAGIC Collaboration has provided new observational data pertaining to the TeV J2032+4130 
gamma-ray source (within the Cygnus OB2 region), for energies _E 7 > 400 GeV. It is then appropri- 
ate to update the impact of these data on gamma-ray production mechanisms in stellar associations. 
We consider two mechanisms of gamma-ray emission, pion production and decay (PION) and photo- 
excitation of high-energy nuclei followed by prompt photo-emission from the daughter nuclei {A*). 
We find that while the data can be accommodated with either scenario, the A* features a spectral 
bump, corresponding to the threshold for exciting the Giant Dipole Resonance, which can serve to 
discriminate between them. We comment on neutrino emission and detection from the region if the 
PION and/or A* processes are operative. We also touch on the implications for this analysis of 
future Fermi and Cerenkov Telescope Array data. 



PACS numbers: 07.85.-m, 98.70.Sa 

Two well-known mechanisms for generating TeV 7-rays 
in astrophysical sources are the purely electromagnetic 
(EM) one proceeding via synchrotron emission and in- 
verse Compton scattering, and the hadronic (PION) one 
in which 7-rays originate from ir° production and de- 
cay 0, 0. Recently, we highlighted a third dynamic 
which leads to TeV 7-rays: photo-excitation of high- 
energy nuclei, followed by prompt photo-emission from 
the excited daughter nuclei [1, 0j. In this chain reac- 
tion, the nuclei act in analogy to Einstein's relativistic 
moving mirror to "double-boost" eV starlight to TeV en- 
ergies for a Lorentz boost factor > 10 6 . The important 
role played by the Giant Dipole Resonance (GDR) in 
the photo-disintegration effectively suppresses the con- 
tribution to the 7-ray spectrum below a TeV [5[. This 
process (which we have denoted by A*) has been pro- 
posed 0] as a candidate explanation of the unidentified 
HEGRA source TeV J2032+4130 at the edge of the 
Cygnus OB2 (Cyg OB2) association (see Ref. Q for the 
discussion regarding the HESS source HESS J 1023-575 
at the edge of Westcrlund 2 0). This stellar association 
has been known to harbor a large population of massive 
and early type stars [lj| that can provide the required 
UV target density. 

Recently, the MAGIC Collaboration has reported new 
TeV 7-ray data from this region [llj. Thus, it is of in- 
terest to expand previous analyses to include this new 
data. In this paper, we present a unified analysis of the 
A* and PION mechanisms for gamma-ray production in 
Cyg OB2. This combined analysis indicates the relative 
importance of these two mechanisms. Extrapolation to 
lower energies then allows us to make predictions within 
reach of the Fermi mission [l2l |. 

TeV gamma-ray data from this region has also been 
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reported by the Milagro Collaboration However, 
inclusion of these data in our study will require further- 
analysis to distinguish contributions from the HEGRA 
source and diffuse interactions expected from the larger 
region of observation. Thus we postpone the considera- 
tion of the Milagro data. 

The critical parameters for the PION and A* mech- 
anisms are the ambient hydrogen gas density and the 
UV photon background, respectively. These parameters 
at present are subject to considerable uncertainty. It 
is a further goal of this work to see whether the present 
and evolving gamma-ray data can meaningfully constrain 
these parameters. 

The stellar distribution of Cyg OB2 reveals a rather 
regular and almost circular density profile with the 
center located at (a, S) = (20 h 33 m 10 s , +41°12') and 
with a pronounced maximum slightly offset at (a, S) = 
(20''33" l 10 s ,+41 o 15.7') Star counts show that 50% 
of the members are located within a radius of 21', and 
90% within a radius of 45' around the center. By in- 
tegrating the radial density profile, after subtraction of 
the field star density, the total number of OB stars is 
found to be 2600 ± 400, with a O-type star population 
of 120 ± 20. This suggests that the total mass of the as- 
sociation is about 1O 4 M . Distance determinations set 
the proximity of the Cyg OB2 to d ~ 1.7 kpc [3]. At 
such distance, the inner 21', with half the total number 
of objects, results in a physical radius of i?; n ~ 10 pc, 
with i? ut ~ 30 pc being the radius of the association. 
Projected onto the sky at the distance of Cyg OB2, the 
HEGRA/MAGIC signal from TeV J2032+4130 @, 03 
was observed only in a 3 pc radius cell at the edge of the 
association. With the same angular radius, there are a 
total of ~ 14 cells in the core of the association. The flux 
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in each of these cells is bounded ~ 1% of Crab, about 3 
times less than that of the TeV J2032+4130 cell. On the 
other hand, th e ag e of the stellar association is supposed 
to be 2-4 Myr [jjf . 

As mentioned above, the prediction of the PION 7-ray 
yield is subject to uncertainty in the ambient gas density 
(as well as the ambient cosmic ray flux). We will find 
that for gas densities larger than 0.1 cm -3 in the vicinity 
of the source, the PION mechanism will dominate, with 
the A* mechanism assuming dominance for smaller den- 
sities. A considerably higher density, njj = 30 cm -3 , has 
been estimated from observations of the CO J = 1 — > 
rotational transition [l9j . If taken as representing an 
average over the core of the association, this value im- 
plies a hydrogen gas mass of 3000 Mq, which is about 
30% of the mass found in stars. Arguments have been 
given [2^ that this estimated density [l!| should be in- 
terpreted as an upper bound. Among concerns are the 
CO-H2 conversion factor [2l|, and the size of the region 
used for averaging. More recently, a thorough analysis of 
the region, with higher angular resolution observations of 
12 CO and 13 CO has been presented (22j. Of interest in 
these results is a significant 13 CO cavity right at the TeV 
source position. This is reminiscent of a formerly found 
IR void. Only one 13 CO clump is seen at a position con- 
sistent with the projected position of HEGRA source. 
No claim is made as to any physical connection between 
the clump and the cavity. Although the clump is massive 
(claimed mass at 337 Mq), its size is significantly smaller 
than the size of the extended TeV source. If gamma-rays 
were to originate from interactions with just this clump, 
there is no reason for the source to appear extended in 
instruments such as MAGIC [H] or VERITAS 0. 

We adopt as our fiducial density the low estimate «h ~ 
0.1 cm -3 , which agrees with the analysis in Ref. [2p| . 
This choice allows us to illustrate the crossover point of 
PION dominance versus A* dominance. It is conservative 
and wise at this point to depend on future experiments 
(discussed below) to resolve the issue of PION versus A* 
dominance. 

First we discuss the PION mechanism for 
TeV gamma-ray production. The emissivity (num- 
ber/volume/time/energy) of neutral pions resulting 
from an isotropic distribution of highly rclativis- 
tic nuclei having a power-law energy spectrum 
dn A {E N )/dE N = N A (E N /E )- a , colliding with 
ambient hydrogen, is given by [l[ 
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where Na is the normalization constant with units 1 / vol- 
ume/energy, Eo is set to 1 TeV, E^{E n o) is the min- 
imum energy per nucleon required to produced a pion 
with energy E„o, and da A {E^a, Epf)/dE^o is the differ- 
ential cross section for the production of a pion with en- 



ergy E„o in the lab frame due a nucleus A of energy per 
nucleon E^ = Ea/A colliding with a hydrogen atom at 
rest. The differential cross-section can be parametrized 
by 

(E„o ,E N )~^- xF n0 (x,E N ), (2) 

where x = E„o/Epf, ctq = A 3 / 4 <tq provides a scaling 
of the cross-section with the atomic number [25] , Co — 
(34.3 + 1.88L + 0.25L 2 ) mb, and F„o (x, E N ) = dN„o jdx 
is a fragmentation function. We take 
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with 13^ = a+0.25, (3 = 0.98/Va, r = 2.6/aA a = 3.67+ 
0.83L + 0.075L 2 , and L = \n(E N /TcV) |2f|. Because 
isotropy is implied injl]), it is straightforward to obtain 
the 7-ray emissivity [6j ; it is 
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where E^ n (E 7 ) = E 1 + ml/{AE^). 

Before proceeding to the A* mechanism, we pause to 
compare Eq. ([3]), which is a functional fit to the out- 
come of numerical simulations obtained with the SIBYLL 
event generator [2711 , to data collected at Tevatron by the 
CDF detector [28[. The results of simulations leading to 
Eq. have been reported at nucleon energies of 0.1 TeV 
and 1000 TeV. The latter corresponds to a center-of-mass 
energy of y/s ~ 1410 GeV. The CDF group at Tevatron 
has measured the charged pion spectrum for pseudora- 
pidity \rj\ < 3.5, at cm energies of 630 GeV and 1800 GeV. 
They have provided a fit over the energy range of inter- 
est, quadratic in ln[s/GeV 2 ], with \ 2 = 0.72 for three 
degrees of freedom 



dN, 
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- = (0.023 ±0.008) In 2 s- (0.25 ±0.19) Ins 
+ (2.5 ±1.0), (5) 



valid for 77 = 0. Taking into account that F v + ~ F n - ~ 
F^o , that the spectral dependence on 77 is mild, and that 
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we find that the results in Eq. ([3]) agree remarkably, 
within one standard deviation, to the CDF fit to their 
Tevatron data. 

Now we discuss the A* mechanism of TeV gamma- 
ray production. The photo-excitation (or photo- 
disintegration) rate for a highly relativistic nucleus with 



energy E = jArriN (where 7 is the Lorentz factor) prop- 
agating through an isotropic photon background with en- 
ergy e and number-density spectrum n(e) is [2!| 
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where cr^(e') is the cross section for photo-disintegration 
of a nucleus of mass A by a photon of energy e' in the 
rest frame of the nucleus. 

We assume that n(e) results from thermal emission 
of the stars in the whole Cyg OB2 association, out to 
-Rout ~ 30 pc. We model the association with half 
of the stars uniformly distributed in the inner region, 
i?i n ~ 10 pc, and the other half uniformly distributed in 
the outer shell, i.e. the density of stars in the inner re- 
gion is (i?out/-Rin) 3 — 1 ~ 26 times that in the outer shell. 
To reproduce the size and position of the source of the 
HEGRA signal, the photo-disintegration must occur in a 
region of radius r ~ 3 pc at the edge of the inner part 
of the association, R < R m . The average photon density 
in this region reflects both the temperatures To and Tb 
of the O and B stars, respectively, and dilution resulting 
from inverse square law considerations. The resulting 
photon density is 
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where N ( B j is the number of O (B) stars , i?o(B) ls the 
O (B) star average radius, and 
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is the Bosc-Einstcin distribution of photons emitted from 
a star at temperature Tq(b)- The factor 47/4 is a con- 
sequence of averaging the inverse square distance within 
this distribution for the density and the region where the 
reaction takes place jf|. It is clear, however, that within 
the 3 pc HEGRA hot spot the concentration of stars 
would be above average, and thus hereafter we take as a 
fiducial value for n HEGRA (e) = 1.7n*(e). The 1.7 factor 
encapsulates an uncertainty of ~ 1 to 2.5 (3(j. The re- 
sulting photo-disintegration rate Ra for the value of this 
density will be denoted by R^ EGRA . 

In Fig. [T] we show the dependence on the Lorentz fac- 
tor of i?™ GRA i f° r the stellar ambiance described above. 
For the O stars we have taken No = 130, a surface tem- 
perature To = 40000 K, and radius Ro = 19 R®; for 
the cooler B stars we assign T B = 18000 K, N B = 2470, 
and radius Rb = 8 R@. The numbers No and N B are 
consistent with the Cyg OB2 data discussed in the intro- 
duction. The cross section has been calculated in both 
the Narrow- Width Approximation (NWA) and the more 
accurate dipole approximation, for the nuclear parame- 
ters given in Ref. Q. In the calculation that follows we 
adopt the more accurate dipole form for the cross section. 

The low-energy cutoff on i?^ EGRA is evident in Fig. [1] 
This cutoff will be mirrored in the resulting gamma-ray 




FIG. 1: Photo-disintegration rates of 56 Fe, 28 Si, and 4 He for 
the HEGRA hot spot. We have approximated the cross sec- 
tion of the GDR by a dipole (solid-lines) and by a single pole 
of the NWA (dash- lines). 



distribution. Notice that the NWA, which we do not 
use below, overestimates the severity of the low-energy 
cutoff. 

The energy behavior for photons in the 0.5 — 10 TeV 
region of the HEGRA and MAGIC data is a complex con- 
volution of the energy distributions of the various nuclei 
participating in the photo-disintegration, with the rate 
factors appropriate to the eV photon density for the var- 
ious stellar populations. Approximating the 7-ray spec- 
trum as being monochromatic with energy E' A (in the 
nucleus rest frame), the emissivity becomes [(| 
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where E 1 is the energy of the emitted 7-ray in the lab, 
and Na, which we take to be 2 [3l| . is the mean 7-ray 
multiplicity for a nucleus with atomic number A. 

It is important to note that the same nucleus source 
density dnA/dE^ is present in the A* emissivity (|10[) and 
in the PION emissivity (J3J) (via ((T|)). Thus, a comparison 
of the two mechanisms will depend only weakly on the 
exact features of dnA/ 'dE^ ■ 

The differential photon flux at the observer's site (as- 
suming there is no absorption) receives contributions 
from both mechanisms, PION and A*. The result is re- 
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FIG. 2: Eyeball fits to HEGRA and MAGIC 7-ray spectrum. 
We have assumed an iron nuclei population with a = 2 and 
E' Fe = 2 MeV. The thick straight red line is a representa- 
tive fit to the combined spectra assuming the PION process 
only with n H = 2 cm" 3 and N Fc = 5 x 1CT 12 cm" 3 TeV" 1 . 
The thick solid blue curve is a similar fit combining both 
the A* mechanism (thin solid blue line) and PION mech- 
anism (dashed blue line) for tih = 0.1 cm" 3 and A/fc 
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servations, a squared-energy weighted differential flux of 
dR 



E: 



dE^ 



3.07 x 1CT "crrT 2 s" 1 TeV (13) 



is obtained. This is nearly two orders of magnitude 
above extrapolation of the HEGRA/MAG1C measure- 
ment at 100 GcV. The presence of low energy power- 
ful sources can clearly dominate the flux at lower en- 
ergies [331 ] . Thus, the normalization inferred from this 
low-energy data can grossly over-estimate the predicted 
flux at 100 GeV (see Ref. [34[ for a thorough discus- 
sion of these issues). A prime candidate source for the 
low energy radiation is a pulsar (with spin-down power 
2.6 x 10 35 crgs _1 ), which coincides (within errors) with 
the position of TeV J2032+4130 (4' displacement) [35[. 
The observed GeV emission can plausibly be ascribed to 
electron acceleration in the magnctosphcrc. This radia- 
tion is exponentially cut off in the TeV region. However, 
such very high energy radiation, which is the focus of 
the present paper, could possibly be associated with in- 
verse Compton scattering of the electrons which power 
the pulsar wind nebula, if such exists (see as an exam- 
ple the case of HESS J1825-137 [36|, and many others in 
the recent literature). Alternatively, the TeV radiation 
can originate in the OB association via the A* and PION 
mechanisms. 



lated to to the 7-ray emissivity as 

^(E-y) = ^ [Qf(25 7 ) + Qf(E y )] , (11) 

where Vdi s is the volume of the source region and d 
is the distance to the observer. In Fig. [5] we provide 
a some eyeball fits (thick solid lines) to the combined 
HEGRA/MAGIC 7-ray spectrum, obtained from inte- 
grations implicit in the two emissivities in Eq. (fTT|) . 

The fits are for an iron nuclei population with spec- 
tral index a = 2 and for an average energy of the pho- 
ton (in the nuclear rest frame) emitted during photo- 
emission E' lFe = 2 MeV [sj. The solid thick blue 
curve is a fit using both the A* mechanism (solid blue 
thin line) and PION mechanism (dash blue thin line), 
with njj equal to our fiducial value, 0.1 cm -3 , and 
N Fc = 3 x 10" 11 cm" 3 TeV -1 . The red thick straight 
line is a representative fit to the combined spectral data 
assuming the PION process only with nn = 2 cm -3 and 
iVpo = 5 x 10" 12 cm -3 TeV -1 . For the iron nuclei popula- 
tion assumed in the fits, the target gas density for PION 
dominance at all energies is njj ^0.5 cm" 3 . 

Additional data is becoming available from observa- 
tions of the Fermi satellite. A preliminary measurement 
in the Cygnus region yields an integrated gamma-ray 

flux m 

F 7 (l - 100 GcV) ~ 3.07 x 10~ 8 cm" 2 s" 1 . (12) 
If one naively assumes a spectrum oc E~ 2 for their ob- 



We now conclude with a discussion of our results: 

• From Fig. [5] it is apparent that the combined 
HEGRA/MAGIC data can be fit with only the 
PION mechanism in operation. Such a fit applies 
if the gas density nn is larger than 2 cm" 3 . For 
0.05 cm" 3 < uh ;$ 2 cm" 3 , a combination of PION 
and A* can provide a satisfactory fit to the data, 
whereas for nn < 0.05 cm" 3 a good fit to all the 
data can be obtained using only the A* mechanism. 

• The average energy of the photon (in the nucleus 
rest frame) emitted during photo-emission has been 
taken as 2 MeV. This is appropriate for iron nuclei. 
If nn < 0.05 cm -3 , a better fit to all the data can 
be obtained using only the A* mechanism with a 
lower average energy of 1.5 MeV. 

• For low gas densities, the spectral features charac- 
teristic of the A* mechanism become visible. These 
are best described as a broad bump in the spectrum 
in the region 1 — 10 TeV. 

• In completing the explanation of the HEGRA and 
MAGIC signal, there remains one issue to address - 
the signal was observed only in a 3 pc radius cell at 
the edge of the inner association. At this point in 
our understanding we can provide only qualitative 
remarks. One possibility is an increased density 
of very hot OB stars in the TeV J2032+4130 cell, 
which provide efficient trapping and accelerating 
conditions for the nuclei, as well as a hot photon 
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background. Indeed, a recent estimate [37J indi- 
cates around 10 O stars in the region of the source, 
a number which is a factor of 3 larger than that 
expected on the basis of a uniform population. 

• If the energy spectrum of cosmic electrons oc E~ 2 
(with an exponential cutoff at 40 TeV), the data 
can also be explained by inverse Compton scat- 
tering of these electrons on the cosmic microwave 
background photons [ll[ ■ The EM explanation can 
only accommodate the data if the Compton peak is 
matched to the energy range of HEGRA/MAGIC 
detection, a possibility allowed within errors. 

• We expect a flux of TeV v^, from both the A* 
(via neutron decay followed by oscillations and 
PION (via 7T^~ decay) [2(j mechanisms. Allowing 
about one muon neutrino per photon after oscilla- 
tion, we expect about 1.2 events/yr at IceCube with 
a background from atmospheric neutrinos of about 
1 event/yr (3{|. However, it is possible that this 
event rate can be considerably enhanced by emis- 
sion from the additional 3-pc cells in the association 
(which will not be resolved by future neutrino de- 
tectors). The signal enhancement can amount to as 
much as a factor of about 5 due to the emission at 
the upper limit value set by gamma-ray observation 
from each cell in the rest of region (e.g., MILAGRO 
measurement in a region centered in the HEGRA 
region but ten times larger [l3[ and MAGIC up- 
per limit in the direction of Cyg X3 [ll[ , which ap- 
proximately coincides with that of Cyg OB2). Such 
accumulation could make the source visible in neu- 
trinos at IceCube. We also note that absorption 
of gamma-rays at the center of the association (see 
e.g., Ref. [401 ) could be relevant, implying an even 
higher neutrino flux from some cells. Observation 
of a neutrino flux from the HEGRA/MAGIC source 
could disqualify an EM explanation of the origin of 
the gamma-rays, at least for this source. 



• The future Cerenkov Telescope Array (4l| will pro- 
vide stronger spectral discrimination between the 
PION and PION+A* mechanisms. This telescope 
is projected to have a factor > 10 larger sensitivity 
than MAGIC/ VERITAS at TeV energies. It will 
also cover the lower GeV energy region (down to 
tenths of a GEV) where the A* mechanism is sup- 
pressed, thus allowing the possibility of comparing 
the two mechanisms with a single data set covering 
the entire energy region of interest. Due to its su- 
perb angular resolution (expected perhaps at a fac- 
tor of 2 or 3 better than that of MAGIC) and field 
of view (several degrees), it will become the ideal 
instrument to distinguish emission components in 
this energy region, and to study morphology of the 
radiation from TeV J2032+4130. 
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